Fluctuation spectrum of fluid membranes coupled to an elastic meshwork: 
jump of the effective surface tension at the mesh size 
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We identify a class of composite membranes: fluid bilayers coupled to an elastic meshwork, that 
are such that the meshwork's energy is a function F e i[^] not of the real microscopic membrane 
area A, but of a smoothed membrane's area A$, which corresponds to the area of the membrane 
coarse-grained at the mesh size £. We show that the meshwork modifies the membrane tension 
a both below and above the scale f, inducing a tension-jump Act = dF c \/dA^. The predictions 
of our model account for the fluctuation spectrum of red blood cells membranes coupled to their 
cytoskcleton. Our results indicate that the cytoskeleton might be under extensional stress, which 
would provide a means to regulate available membrane area. We also predict an observable tension 
jump for membranes decorated with polymer "brushes" . 

PACS numbers: 87.16.Dg, 68.03.Cd, 87.16. Gj, 68.47.Pe 



Random surfaces [l| play an important role in many 
areas of physics, spanning from biophysics [2j and chem- 
ical physics to high-energy physics Q. In the last 
decades, there has been considerable interest in the 
properties of fluid bilayer membranes composed of am- 
phiphilic molecules in water. It is well-known [f| that 
their fluctuations are dominated by the membrane bend- 
ing rigidity, k, which leads to a spectrum proportional 
to (k<j 4 + aq 2 )- 1 HUH. Here, a is an effective sur- 
face tension, i.e., an adjustable thermodynamic parame- 
ter arising from the constraint of constant surface area, 
which is usually several orders of magnitude smaller than 
the surface tensions of ordinary liquid interfaces 0, ^| . 
Additional complexity arises if the membrane interacts 
with other systems (e.g., a rigid substrate, another mem- 
brane, a network of polymers or filamen ts) fill H2I fl3 | . 
or contains inclusions (passive or active) [lj, [l5j . 

In this paper, we identify a class of composite mem- 
branes, i.e., of membranes coupled to an external sys- 
tem (like the cytoskeleton of red blood cells (lit), for 
which the coupling energy can be described, in a first ap- 
proximation, by a function of the membrane area coarse- 
grained at a characteristic length scale £ (the mesh size in 
the cytoskeleton case). For such systems, we show that 
the effective membrane tension should exhibit a jump of 
amplitude Act at the scale £. In other words, the fluc- 
tuation spectrum is proportional to (k<? 4 + <r < q 2 ]~ 1 for 
q < and to («g 4 + a > q 2 )~ 1 for q > As we shall 
see, Act = a < — a > is directly related to the coupling 
between the membrane and the external system. Un- 
derstanding the scale dependence of the effective surface 
tension, and in particular its value at short length-scales, 
is important in a number of phenomena, e.g., membrane 
adhesion [fR fl8j |. cell fusion [lfl Efil ] and other micro- 



scopic biological mechanisms, such as endocytosis |21l . 
It should also help interpret experiments on composite 
membranes [22, yji |2J] , which attempt to determine the 
value of k by fitting the fluctuation spectrum. 

We shall consider two different systems for which our 
model predicts a jump in surface tension: (i) membranes 
decorated with polymer "brushes" |25j , and (ii) red blood 
cell membranes interacting with their attached cytoskele- 
ton |l6|. Polymer decorated membranes are useful in 
providin g a sterical stabilization for liposome drug car- 
riers [2ri l27j: we expect that better understanding their 
membrane tension should yield new insight in their sta- 
bility and mechanical properties. Here, we shall mainly 
focus, however, on the composite membrane of the red 
blood cell (or erythrocyte). The cytoskeleton of red blood 
cells is a roughly triangular flexible protein meshwork, 
composed of spectrins attached to the bilayer by associ- 
ation with integral membrane proteins [16j. The result- 
ing mechanical strength enables red blood cells to un- 
dergo large extensional deformations as they are forced 
through narrow capillaries. In the last two decades, the 
red blood cell membrane has been thoroughly studied 
and a rich des crip tion of its attached cytoskeleton has 
been obtained [23, H^, |3(| • Very recently, the spectrum 
of the short wavelength membrane fluctuations of red 
blood cells [22| has been successfully fitted by includ- 
ing the effect of confinement from the cytoskeleton |24| : 
however, this approach leads to an unexplained, rather 
abrupt change in the value of the membrane tension. We 
shall show that our model explains quantitatively this 
change, provided that the cytoskeleton is under (exten- 
sional) stress, which is an important issue currently de- 
bated in cell biology 0, ■ Specific implications con- 
cerning red blood cells as well as possible experiments on 
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polymer decorated membranes will be discussed at the 
end of the paper. 

Let us introduce our model by considering a membrane 
(embedded in the Euclidian three-dimensional space) 
connected to an elastic mcshwork. The meshwork is 
composed of N Hookean springs (spring constant A, free 
length £q) and is attached to the membrane through its 
nodes. The shape of the membrane is specified by giving 
the position R(s x , s 2 ) of a membrane element as a func- 
tion of two intrinsic coordinates. The total free energy of 
the system may be written as F[R] — F mem b [R] + F e i [R] , 
where F mem (,[R] is the membrane free energy and 
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is the meshwork elastic energy. Since the meshwork can- 
not follow the wiggles of the membranes of wavelength 
shorter than the mesh size (see Fig. 0, the spring elon- 
gation £ [R] actually depends on the membrane conforma- 
tion. To simplify, we have assumed that £[R] is identical 
for all the springs, which is reasonable since the mem- 
brane is fluid. The mesh size, i.e., the spring length, can 
thus be expressed as a 
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where A^[R] is the area of the membrane coarse-grained 
at the mesh size, and g is a constant depending on the 
meshwork topology (e.g., g = 2^3 for triangular mesh- 
work). Determining £[R] is the two-dimensional analo- 
gous of measuring the length of a fractal coast (e.g., the 
coast of Brittany |3j|) with a straight ruler, and deter- 
mining the ruler length that allows to go through the 
whole coast in N steps. 

Solving Eq. (J2J for £ in the general case of an arbitrary 
membrane shape is a rather difficult task. In order to 
simplify the problem, let us consider the case of a quasi- 
spherical membrane of fixed volume V and fixed area 
A. If V is expressed as V = §7i\Ro, the area can then be 
written in the form A = A+SA where A = 4ttRq ^> 5 A is 
called the projected area. We now consider a very small, 
almost planar portion a = eA of the projected area, that 
still contains a rather large number n — eN of springs. 
We represent the local membrane shape by its fluctuation 
modes it q above the local projected area. Locally, Eq. 
can be expressed as 



-£ 2 M = at[u] = a + 5 
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where rs eA^ is the local coarse-grained area. Consis- 
tently with the second-order expansion, the upper limit 
£ _1 [it] in the above sum can be replaced by its zeroth 
order approximation: 
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FIG. 1: Section of a fluctuating vesicle (thick line) to which is 
attached a meshwork of springs (dashed lines), with common 
length £. The actual membrane area is A and the mesh area is 
A(_, which corresponds to the membrane area coarse-grained 
at f. 



The elastic energy / e i = |nA(£[u] — £q) 2 » eF c \ of the 
small portion under consideration is then given by 
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To proceed further, we introduce the standard form of 
the local membrane elastic energy (/mpmb ~ efmemb)) 
which is given, to quadratic order, by |3j 



(7) 



Here 7 is a confining potential, a an effective tension yet 
to be determined (a Legendre transform multiplier con- 
strainingthe total membrane area), and k is the bending 
rigidity 5]. From the total energy f c i[u] + / mC mbM, we 
deduce, using equipartition, the fluctuation spectrum: 
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The adjustable parameter a is determined from the re- 
quirement 5? 2 (l' lt q| 2 ) — e8A. It follows that the de- 
composition of the low wavevector tension, cr < , and of 
the high wavevector tension, cr > , in terms of a and Act 
(ct < = ct + Act and ct > = a) is quite arbitrary: only 
the values of ct < and ct > really matter. Note that the 
position ^q 1 of the jump only enters as a rough char- 
acteristic scale in our model; a more refined calculation 
would be required to determine the wavevector width of 
the jump. We have also omitted the corrections to the 
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bending rigidity originating from the meshwork, which 
may also be scale-dependent. 

The results obtained above are based on the fact that 
the meshwork energy is actually a function of -A^R] 
through Eqs. QJ and (J2J. The same conclusions con- 
cerning the effective tension therefore hold for any com- 
posite membrane for which the coupling may be ap- 
proximated by a function F C \(A^) of the membrane area 
coarse-grained at a characteristic scale £. The general 
expression giving the tension jump occurring at £ o ~ £ is 
simply 
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Before making predictions concerning other systems, 
such as polymer decorated membranes, let us discuss in 
more details the cytoskeleton case. 

Inspired by the recent work of Gov et al. [24| , we have 
re-analyzed the data of Zilker et al. on the fluctua- 
tion spectrum of normal red blood cells in the wavelength 
range 300 nm < A < 6000 nm. By fitting these data, the 
authors of Ref. [24| have demonstrated the necessity to 
introduce a confining potential, arising from the presence 
of the cytoskeleton. They acknowledge, however, that the 
data at very short scales are better fitted by a value of the 
membrane tension much smaller than the one obtained 
at large scales. Since the cytoskeleton is well described 
by a meshwork of elastic springs [lil l32| , we believe that 
the theory put forward in this paper provides an explana- 
tion to this jump in surface tension. Fluctuation spectra 
are usually fitted via an effective wavevector-dependent 
bending rigidity K q defined by 

<si 2 > = (io) 

According to Eq. JSJ), this quantity should roughly ex- 
hibit a discontinuity at the mesh wavevector qo = £o ■ 
This is indeed the evidenced by the solid line 

fit in Fig. |3 The position of the discontinuity being 
fairly marked at ^q 1 ~ 125 nm -1 , and the value k ~ 
2 x 10~ 20 Nm being already well documented [10113, 
we have fitted the data to Eq. JHJ for the parameters 7, 
ct < and <j > only, obtaining 7 = (5.2 ± 0.6) x 10 6 Jm~ 4 , 
ct< = (1.5 ± 0.2) x 10~ 6 Jm~ 2 and ct> = (-0.8 ± 
1.2) x 10~ 7 Jm~ 2 . According to our model, Act = 
ct < — a > ~ 1.6 x 10~ 6 J m~ 2 should be given by Eq. 10. 
Indeed, assuming £0 > £q (and yet £ ~ £0) we deduce 
A w 10~ 6 Jm~ 2 , which agrees well with the measured 
cytoskeleton elastic constants j^ij. From the positive 
sign of Act, our analysis predicts that the cytoskeleton 
is under extensional stress: this important issue would 
require more experiments to test. Finding ct > essentially 
zero, or even negative, is also quite remarkable. It sug- 
gests that one important role of the cytoskeleton might 
be to regulate, via its elastic stress, the membrane ten- 
sion at length-scales shorter than the cytoskeleton mesh, 
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FIG. 2: Fit (solid line) of the measured fluctuation spectrum 
of a normal red blood cell (squares) from the data of Zilker et 
al. |2^| . The dashed line extrapolates the large length-scale fit. 
The discontinuity and the corresponding tension jump occurs 
at a length-scale ~ 125 nm comparable to the cytoskeleton 
mesh. 



in order to produce either large fluctuations, or even a 
short-scale buckling. This new type of area regulation 
(see Ref. [35j) might prove useful to allow for large de- 
formations. 

Another system to which our model can be applied is 
a membrane bearing an anchored polymer brush. For a 
fixed number p of polymers in good solvent conditions, 
we have 1251 



M, 5 / 3 (l-2 X ) 1/3 £ 5/6 
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where N denotes now the index of polymerization, b the 
monomer length, \ t ne Flory parameter, and £ the graft- 
ing density. As illustrated in Fig. [31 the relevant area 
defining the grafting density is not the total area A, but 
the area coarse-grained at the blob scale. Hence, 
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Expressing as in Eq. J5J the brush energy in terms of the 
local coarse-grained area, or equivalently, using Eq. 10, 
we obtain 



Act 



-iVfc B T& 5 / 3 (l-2 X ) 1 / 3 (£) 11/6 . (13) 



This actually corresponds to the brush osmotic pressure 
integrated over the brush height: Act w — II osm H ~ 
-k B T(p/A) 3 / 2 x bN^p/A) 1 / 3 . Note that Act < 0, i.e., 
the brush reduces the large-scale tension with respect to 
the short-scale one. It actually draws some of the mem- 
brane area stored at length-scales shorter than the blob 
size in order to lower its grafting density. Let us esti- 
mate Act in typical experiments [l2j: with N = 100, 
6=1 nm, and a (projected) area per chain of 100 nm 2 , 
we obtain | Act| w 10~ 4 J m~ 2 , which compares with ordi- 
nary membrane tensions. At length-scales larger than the 
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FIG. 3: Sketch of a polymer "brush" anchored on a fluctu- 
ating membrane. The actual grafting density is the number 
of polymers divided by the membrane area (dashed line) 
coarse-grained at the blob size £. 

brush height H , the brush elasticity also yields a renor- 
malization of k |3(|. Note that for distortions wavevec- 
tors in the interval H^ 1 < q < the brush will ac- 
tually heal over ~ q^ 1 , implying that the tension jump 
will probably be spread over this interval Our re- 

sult agrees with a recent direct calculation of Bickel et 
al. |38j . showing that membranes bearing polymers in the 
"mushrooms" regime should possess an excess tension 
~ 2/cbTS at scales shorter than the gyration radius R g . 
Indeed, at the grafting density £ ~ which crosses- 

over to the brush regime, this gives Act ~ — LI osm iJ, with 
LI osm ~ kBT/Rg and H ~ R g . A possible way to test our 
prediction Eq. I|13|l would be to measure, on a given vesi- 
cle decorated with a polymer brush, the variation of the 
large-scale membrane tension induced either by a modi- 
fication of the solvent quality or by a chemical alteration 
of the index of polymerization. 

To conclude, the jump of the effective surface ten- 
sion evidenced in this paper arises from a competition 
between the free energy of the membrane, which re- 
quires a certain amount of excess area to be stored 
in the short wavelength modes, and the free energy of 
the coupled elastic system, which is lowered if either 
less or more excess area (depending on the system) is 
transferred to the large-scale modes. Very generally, a 
wavevector-dependent membrane tension should there- 
fore occur whenever an external system coupled to the 
membrane produces a scale-dependent area transfer. 

We acknowledge fruitful discussions with A. Ajdari, 
F. Gallet, and P.-G. de Gennes. 
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